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Abstract 

We have made a detailed and quantitative study of dilepton production via 
bremsstrahlung of a virtual photon during pion-pion collisions. Most cal- 
culations of electromagnetic radiation from strong interaction processes rely 
on the soft photon approximation (SPA). The conditions underlying this ap- 
proximation are generally violated when dilepton spectra are calculated in 
terms of their invariant mass, so that an approach going beyond the SPA 
becomes necessary. Superseding previous derivations, we derive an exact 
formula for the bremsstrahlung cross section. The resulting formulation is 
compared to various forms based on the SPA, the two-particle phase space 
approximation and RiickFs formula using a relativistic One Boson Exchange 
(OBE) model. Within the OBE approach, we show that approximations to 
the bremsstrahlung dilepton cross sections often differ greatly from the ex- 
act result; discrepancies become greater both with rising temperature and 
with invariant mass. Integrated dilepton production rates are overestimated 
by Riickl-based approximations by factors 1.5-8.0. The largest discrepancies 
occur for the reaction ir + ir + — * ir + ir + £ + l~ , where such approximations over- 
estimate the exact rate by factors ranging from 2 to 30 for invariant masses 
between 10 and 500 MeV. Our findings, combined with recent estimates of 
the Landau- Pomeranchuk effect, indicate that bremsstrahlung dileptons rates 
in ultrarelativistic heavy ion collisions should be even more suppressed than 
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had been thought before. 
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I. INTRODUCTION 



Dileptons produced in heavy ion collisions interact only electromagnetically with their 
hadronic surroundings. Especially the early phases of such collisions, when dilepton produc- 
tion rates are largest, can therefore be explored with lepton pairs. Some of the initial ideas 
on photon [l| and dilepton |2| [|J emission have met with considerable interest, and this 
area has now expanded into a major component of research in ultrarelativistic heavy ion 
collisions [ffl. On the experimental side, dilepton experiments for hadronic collisions ||- || 
have been complemented only recently by nucleus-nucleus experiments at very high energies 
M- fT2|]. Since the spectrum of lepton pairs and of real photons has been proposed as a 



signature of the quark-gluon plasma 0, it is imperative quantitatively to understand the 
emission mechanisms from the confined sector of QCD, which in this sense constitutes the 
background. 

Guided by theoretical estimates, we can roughly divide the lepton pair production cross 
section, expressed in terms of its invariant mass M, into a number of distinct regions. For 
M > mj/^, the Drell-Yan contribution dominates, while pairs arising from the semileptonic 



decay of DD are also important fll3 |. For lepton pair invariant masses around the light 
vector meson masses (p, u>, and 0), two-body |14| and three-body [|15[] reactions dominate. 
The latter processes also extend their influence to the region < M < mj/^, where the 



quark-gluon plasma may become visible for sufficiently high initial temperatures ||16|| . In 
this paper, we concentrate on the regions of low invariant masses, the so-called "soft" limit 
M < 300 MeV. There, production via bremsstrahlung is expected to play an important role 



and many authors have recently concerned themselves with such reactions p77| - ||28|| . 

Since particle production plays an important role in ultrarelativistic heavy ion collisions, 
we consider only events that arise from microscopic meson-meson dynamics. In high energy 
heavy ion experiments, the meson to baryon ratio is such that the absence of baryons in our 
treatment is not expected to constitute a major hurdle to phenomenology. More specifically, 
we restrict our studies to those of pion-pion bremsstrahlung. 

The production of photons or dileptons necessarily involves 3- or 4-particle final states. 
Except for special cases such as cross sections at fixed q P9| , such final states have universally 
been handled under some approximation. It is our goal in this paper to go beyond these by 
presenting an essentially exact formulation for the ttti — > r n r nt r t' cross section (as a function 
of dilepton invariant mass) and, in the process, to test how far these approximations may 
be trusted. In order to achieve these goals, we formulate the pion-pion interaction in terms 
of a relativistic One Boson Exchange model (OBE). 

Note that we do not attempt to compare directly to experimental data because the 
main point at issue is a comparison of theoretical approaches to bremsstrahlung. To do 
justice to small-M data would require spacetime integration over the fireball region as well 
as inclusion of Dalitz decays and the Landau- Pomeranchuk effect, a substantial task which 
will be postponed to future work. 

Our paper is organized as follows: in the next section, we discuss the set of approxi- 
mations known under the collective banner of the Soft Photon Aproximation (SPA). Next, 
we derive exact formulae for the bremsstrahlung generation of lepton pairs. Section |V| is 
devoted to technical issues having to do with phase space; Sections [V] and [VI| then derive 
cross sections based respectively on various SPA variants and the exact formulation. Our 
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OBE model is presented in Section |VII| ; finally, Section [VIII| contains our numerical results 
in the form of cross sections da(s)/dM and rates dN/d 4 x dM and a discussion on the issues 
raised. 



II. THE SOFT PHOTON APPROXIMATION 

The Soft Photon Approximation (SPA) has been widely used in calculating bremsstrah- 
lung dilepton spectra. It is based on the early realization |30| that the cross section for 
production of low-energy real photons is dominated by the corresponding hadronic ampli- 
tude. In practice, this means neglecting the photon momentum in the hadronic matrix 
element as well as all photon emission from vertices and internal lines (see Section 1IB ). 



A. Limitations 



In order for the SPA to be a valid approximation, two conditions [31| must be met: 



1. The photon energy must be much smaller than the energy of any one of the hadrons 
participating in the scattering, 

?o/#<l> (1) 



2. the process of radiating the photon must be separable from the hadronic interac- 
tion process, i.e. the hadronic and electromagnetic time and distance scales must be 
sufficiently different to permit separate treatment. The range of the hadron-hadron 
interaction is given roughly by b = rriy 1 , the inverse mass of the exchange boson, while 
the distance the hadron may propagate with off-shell energy AE = q Q before emitting 
the photon is about Ax = v/AE = \p\/Eq Q . The SPA is valid only when Ax/b ^> 1, 
or 

q <^m Y \p\/E. (2) 



Implicit in these equations is, of course, a specific Lorentz frame with respect to which the 
energies are measured. For the special case where two hadrons of equal mass m collide, 
these conditions can be re-written in their cms frame as 

q*o « v^/2 , (3) 

<?o < m Y \Jl - Am 2 /s. (4) 

In a simple bremsstrahlung experiment, these limits are easily satisfied by selecting only 
photons or dileptons of low energy in the laboratory frame. An early dilepton paper by Riickl, 
for example [fffil , calculated spectra for y/s = 27 and 53 GeV while restricting transverse 
momentum to below 500 MeV. 

In the complex multiparticle systems formed in the course of nucleus-nucleus collisions, 
however, the situation is more complicated: There are many binary collisions, and their 
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respective cms frames do not generally coincide either with one another or with the overall 
nucleus-nucleus cms frame. A photon that is soft in a particular hadron-hadron cms therefore 
does not have to be soft in the laboratory and vice versa. (This fact has to be taken into 
account when addressing issues related to the spectrum of real photons f33j.) 

In a situation where there are many cms frames, it therefore becomes necessary to look at 
relativistically invariant quantities. For lepton pairs, the extra degree of freedom provided by 
a non-vanishing invariant mass is a suitable variable of choice. Looking at invariant masses 
means, however, that go is no longer fixed but must vary over its full kinematic range, which 
for our example of colliding equal-mass hadrons is given by (see Section |V|) 

s - Am 2 + M 2 

Once the upper kinematic limit for q approaches or exceeds the bounds set by either Eq. 
(0) or ([|), the Soft Photon Approximation clearly fails. 

In Figure 1, we show the three functions (§) and @ and (H) for the case m = m n = 140 
MeV, my = ma ~ 500 MeV and dilepton invariant masses M — 10 and 300 MeV. It is 
immediately clear that the assumptions underlying the SPA are not fulfilled even for small 
M: the kinematic range accessible to go is never much smaller than the limits set by the 
SPA. The situation becomes even worse for larger M. 

Closely related to the above is a common misconception that the SPA must necessarily 
be valid for small invariant masses. The confusion stems from the fact that, because 



M 2 = q 2 -q 2 } (6) 

a small value for go necessarily implies small M. However, the reverse is not true: both 
go and q can be large even while adding up to small M. An extreme counterexample is of 
course a hard real photon with zero invariant mass but large energy. One cannot, therefore, 
a priori expect SPA calculations to coincide with full hadronic calculations even when M is 
small. Section |V111] provides examples of such deviations; see especially the ratios between 
approximations and the exact calculation in Figures 12 and 13 which do not approach unity 



for small M. This point has been emphasized before [19| . 

We are therefore led to conclude that the SPA does not generally apply to rates or cross 
sections measured in terms of the dilepton invariant mass. Correspondingly, doubt must be 
cast on the validity of the SPA in calculations pertaining to such measurements. 

That said, one may still ask by how much SPA-based calculations differ from the more 
complicated true situation: is it a matter of a few percent, or orders of magnitude? Full 
calculations such as presented in Section [VI I| below are cumbersome and time-consuming, 
and it is therefore of value to understand quantitatively to what percentage the SPA approx- 
imation can be trusted. It is this quantitative question that we attempt to cast some light on 
by making a model calculation capable both of "true" answers and an SPA approximation. 

Clearly, a single model calculation such as presented below cannot answer this quanti- 
tative question exhaustively. Nevertheless, its results provide at least an indicator of the 
reliability of the SPA: If, within our simple OBE, dileptons can or cannot reasonably be 
described within the SPA, it is possible or even likely that the same or similar conclusions 
would emerge from other models. 
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B. The SPA as series of approximations 



Rather than being just a single step, the soft photon approximation as implemented in 
a practical context can consist of a number of stages. Starting with the most fundamental, 
these are: 

(A) The virtual photon producing the lepton pair is emitted only by external legs of the 
hadronic reaction; emission from internal parts of the hadronic reaction is neglected 
because the radiation from internal lines forms a sub-leading contribution |31[|. 



(B) The dependence of the hadronic matrix element Aih on the photon momentum q is 
neglected. 

(C) The photon momentum q is neglected in the phase space delta functions. 
In addition, we shall be considering in this paper the effect of 

(D) approximating the virtual photon current by the real photon current, (J M = J2{2Pi ± 
q)/(2 Pi -q±M 2 )-> Epf/(Pi ■ q)\ see Eqs. (0) and ©). 

(A) and (B) have formed the basis for most bremsstrahlung dilepton calculations as they are 
hard to avoid. The restriction of phase space (C), on the other hand, can be circumvented 
fairly easily by inserting later a phase space correction factor (Section |V A ) or by avoiding 
(C) altogether using Eqs. (p5|) , fl36|) and Sections [TV] and [VT|. 



The change from virtual to real photon currents in step (D) is not directly dependent 
on (C) and can be implemented either jointly or separately. We shall be investigating both 
alternatives and determining their impact on the overall dilepton rate. In Section [II IB] and 
Section |V1I| and beyond, we shall refrain from using even the basic assumptions (A) and (B) 



and calculate exact cross sections and rates within the same lagrangian. 



III. BREMSSTRAHLUNG CROSS SECTION 

We wish to find the cross section for a dilepton pair £ + £~ emitted in the semi-elastic 
scattering of two pions a and b into 1 and 2, 

vr a + tt 6 -> 7Ti + tt 2 +£+ + £- . (7) 

The lepton pair forms via a virtual photon 7* with four-momentum q^ and mass q 2 = M 2 
emitted from either the central blob or external legs of the hadronic reaction 

TTa + 7T& -> 7Tl + 7T 2 . (8) 

The leading contribution will come from photon emission from external pion legs. First, 
we follow Lichard [^6j in a leading-terms derivation but place more emphasis on the role 
of phase space; the full formalism supplanting the leading-terms derivation is presented in 
Section I111B. 
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A. Bremsstrahlung from external legs 



Four diagrams contribute to the leading-term cross section, one of which is shown in 
Figure 2. In the notation of Figure 2, the total matrix element can be written as products 



of purely hadronic reactions Ai^ with currents Jf and a leptonic part L„ [34 



where 



M e+e - = e [J£M h {p a - q,p b ,Pi,P2) + Jb M h{p a ,Pb ~ q,P\,Pz) , 

+JiM h (p a ,p b ,p 1 + q,p 2 ) + J 2 M h (p a ,p b ,p 1 ,p 2 + q)] , (9) 



g 

L » = jpU{P-)l^{P+), (10) 



and, in slightly condensed notation, the terms of the pionic current are 

J *- 2p a ,-q-M* ' {U) 

Tlt _ Qi,2(2pi,2 + qV 
2 Pi,2 ■ q + M 2 

with Qi the charge of pion % in units of the proton charge. Note that Jf contains the prop- 
agator for the off-shell pion. Making approximation (B), one now neglects the dependence 
on q of the four hadronic matrix elements; for example 

M h {p a - q,Pb,Pi,P2) ~ M h (p a ,p b ,p 1 ,p 2 ) , (13) 

permitting the usual factorization of Aie+e- into a purely hadronic matrix element and a 
photon-plus-lepton part: 

M e +e- ^ eM h (p a ,p b ,pi,p 2 )J^(p a ,Pb,Pi,P2,q)L fl (p + ,P-) , (14) 
where now the virtual-photon current 

contains the sum over the four diagrams. Since approximations (A) and (B) are valid only 
for photons of small q and by implication for small M, the extra terms in numerator and 
denominator are usually neglected (in what we called approximation (D)) to form the real- 
photon current 

p = E Qi— ~ E Qi— ■ ( 16 ) 

i=1^2 Pi 'I i=a,b Pi-1 

Squaring the total matrix element and summing over lepton polarizations (s+,s_) yields, 



using either current p6 
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E \m 



£+£- 



ai Pby Ply P2)\ J^J L^ v 

,327r 2 a 2 r , , (I -J) 2 ' 



\M, 



M 2 



-J- J 



M 2 



where I = p + — p- and the leptonic tensor is 

8na 



J fJ,U 



E W 



M 4 



(17) 



(18) 



At this point, q is still defined only as the sum of lepton momenta (p + + pJ). 

With F = 4[(p a • pb) 2 — m 2 m 2 } 1 ^ 2 the incoming flux, the unpolarized cross section is an 
integral over phase space of four outgoing particles, 



do~hhe+e- 



E \M 



dH,A 



e+e- 



(2tt) 8 F 



dR A = 5 A {p a + p b - pi - p 2 - p+ - P-) dri dr 2 dr + dr_ 



(19) 
(20) 

where we write dr^ = d 3 pi/2Ei for short. This can be reduced to a t/iree-particle phase 
space integral by changing variables and integrating over redundant degrees of freedom in 
the dilepton cms, 



dr+dr- = -i/l — 



An 2 

M 2 q 4 



(21) 



where the positron solid angle fl + must be kept due to the (/ • J) term, and /i is the lepton 
mass. The dependence of (I ■ J) on Q + is integrated out, 



-J 



M 2 



8tt 



2fi 2 



dQ + = T [ 1 + W 2 M- J 



(22) 



so that, using Eqs. |17H22|, we can express the total cross section in terms of a leptonic factor 
K and a "hadronic plus virtual photon" cross section ddhh-y*, 

da hhi+e - = k(h 2 , M 2 ) da hhl * , (23) 



with 



and where 



k = M 2 t(M 2 



da. 



hh-y* 



Ana 



a 
3tt 

dM 2 
Jl 2 



K(p ,M ) da hhr , 
2^ r 



1 + 



M 2 



1 



4/i 2 
M 2 



'-J 2 )\M } 



dR* 



(2tt) 5 F 



is a function of i/iree-particle phase space 



dR 3 = 6 4 (p a + p b -Pi-P2-q) dn dr 2 dr ( 



g ■ 



(24) 



(25) 



(26) 



With the understanding that its derivation is valid under approximations (A) and (B) only, 
Eq. (|23D is exact and makes no assumptions with respect to phase space. In particular, we 
note that q is still contained in the delta function constraining the phase space in Eq. (f2B|). 
Next, we shall show that Eq. (|23|) is valid even when approximations (A) and (B) are not 
made but that da^* takes on a form different from Eq. (|25|) . 
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B. Emission from all diagrams 



When calculating dilepton cross sections for all possible diagrams, including internal 
emission of photons, it is no longer possible to factorize the matrix element into a hadronic 
part and a current as in Eq. (P). Instead, we consider, for a given pion-pion reaction, the 
full matrix element 

M ( nn ^ nni+r) = M » Lfi = Y / M^, (27) 

m 

where m runs over all contributing diagrams. The sum of diagrams is gauge invariant, 

%W = £ q,M» m = , (28) 

m 

but not necessarily the individual terms. Unlike in the previous sections Ai^ contains, 
besides the hadronic interaction, the photon vertex and where applicable the external pion 
propagator, both of which were previously part of the current. 

Starting again with the cross section in four-phase space (0) and transforming from 
dr + dr_ to dr q dM 2 dfl + using Eq. (pTD, the unpolarized dilepton cross section is 



/ 4u 2 dM 2 dT 1 dT 2 dT g 5(- ■ ■) f ^ . , . l2 
iahh „. = - JL 8(2>r)8 ; 1 J dU + £ |AW 1 . (29) 

Squaring and summing over lepton polarizations, we obtain in analogy to Eq. fll"7|) 

E |a^/ = E^^:%- (30) 

s+s- ™,n 

where as before L au = 87ra(q n q u — l a l u — M 2 g^ v )/M A . In the dilepton cms, l = and the 
second term in L^ u integrates to 

J dQ+l^M^M™ = 4 J dn + ( P+ -M m )( P+ -M* n ) , 

Att 

= -(M 2 -A^)M m -M* n . (31) 

Since q = in this frame, = qa-M-m/qo and hence in the dilepton cms 

Mm'M* n = q» 2 (q»M»J(q„M* n n - M m -M* n , (32) 
so, using gauge invariance fl28|) , 

E M m -M* n (cms) = - E M m -M* n . (33) 

van ran 

Hence 

f*i> e W = (i + 1!) (-E^;) , (34, 
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and the full cross section becomes 



d<J hM +e- = K(M 2 )da hhr , (35) 
which is identical to Eq. (0), but with the virtual photon cross section reading 

dM 2 



da, 



hh-y* 



M 2 



':^ m - m ')wtf- (36) 



Comparing to eq. (p5l), we have the correspondence 4na(— J M J M ) |A4 



y M* 



njx ) 



Equations (^) and (|36|) , together with use of the 3-phase space variables of Section 
|TVl provide a fully covariant formalism capable of exact treatment of both hadronic and 
electromagnetic sectors of a given model. Within the physical limitations implicit in a 
given model, it makes no further assumptions or approximations. In particular, none of the 
approximations listed in Section [11 B| are needed. 



C. Dilepton rates 

In the context of nuclear collisions, one is interested more in overall dilepton rates than 
in fixed-s cross sections. Not attempting to account fully for either quantum statistical or 
flow effects here, we shall make the simplest assumption of a locally thermal Boltzmann gas 
of pions at a (local) temperature T. The rate of pair production of invariant mass M per 
unit of four-volume then is || 

ds\{s,m a ,m b ) TTTo ' ( 37 ) 



d^xdM 2 327r 4 7 v ' a ' bJ (V^/T) dM 2 
where 

Kx,y,z) = (x - y - z) 2 ~ Ayz (38) 



is the basic three-particle kinematic function |35| , which for equal masses (m a = = mi = 
m 2 ) has the more familiar form A(s, m 2 , m 2 ) = s(s — 4m 2 ), and g a b = {2S a + l)(2Sb + 1) is 
the spin degeneracy factor which for pions is unity. K\ is the modified Bessel function. 

IV. THREE-PARTICLE PHASE SPACE 

To complete our preparatory work, we summarize the treatment of three-particle phase 
space in terms of relativistic invariants. The relations of this section will be useful first in 
3-phase space integration of the currents fll5]) and ([16]) and later in the exact treatment of 
bremsstrahlung dilepton production. 

Three-particle phase space has been studied extensively and we here merely outline the 



procedure; for details see Ref. |35[. As shown in Figure 3, the kinematics for a reaction 

a + b — >l + 2 + 3 (39) 
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can be described in terms of the five invariants s = (p a + Pb) 2 , t\ = (pi — Pa) 2 , «i = 
(Pi + P2) 2 , $2 = ij>2 + P3) 2 , and t 2 = (pb — Pz) 2 - For our purposes, we shall be using p 3 and 
q interchangeably as referring to the virtual photon, i.e. E 3 = q and q 2 = m\ = M 2 . The 
other momenta refer, as usual, to the incoming and outgoing hadrons. Inverting, one gets 



2Pa- 


Pb 


= s - m 2 a - m\ , 


2Pa- 


Pi 


= ml + m\ — t± , 


2Pa- 


P2 


= si + tx - t 2 - m\ 


2p b - 


Pi 


= 3-32+^-1711, 


2 Pb - 


P2 


= s 2 + t 2 -t 1 -M 2 


2pi ■ 


P2 


— si — m\ — m\ , 


2 Pa 


■ 1 


= s - si + t 2 - ml , 


2p b 


■ q 


= m 2 b + M 2 - t 2 , 




■ 1 


= s - Si - s 2 + ml , 


2p 2 


■ 1 


= s 2 - mj - M 2 . 



(40) 



The corresponding phase space integral is 

7T 



dR 3 (s) = 



4A x / 2 (s, ml, ml 



dti ds2 ds\ dt 2 



(41) 



where the weighting is given in terms of the Cayley determinant 



B 



1 

1 

1 m\ 



mt. 



s 2 ti m\ 



I M l t 2 si 
1 s 2 M 2 ml s 

,2 „ „ ^2 



1 ti t 2 mi m 2 



m" a 



(42) 



1 m\ S\ 

The Cayley determinant is quadratic in any of its arguments; specifically, we use the form 



B = \{s,s 2 ,m\) (4 -t 2 )(t 2 -q), 
where the kinematic limits on t 2 are given by 



(43) 



m 2 b + M 2 



X(s, s 2 ,m\) 



2s 

s+s 2 -m\ 
,2 1 ™2 



5+52—^1 S — Si+M 2 

2s 2 s 2 -m 2 +M 2 




s-ml+m'l s 2 -ti+ml 
G(s, ti, s 2 , m a , m b ,m 1 )G(s 1 , s 2 , s, m 2 , m x , M 



1/2 



(44) 



\(s,s 2 ,m\) 

written in terms of yet another determinant and the basic four-particle kinematic function 



G{x,y,z,u,v,w) 



1111 

1 v x z 
1 v u y 
1 x u w 
1 z y w 



(45) 
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The kinematic limits on s± are 

sf = s + M 2 - — (s + s 2 - m\)(s 2 + M 2 — m 2 ,) 
2s 2 

± t^-A 1/2 (s, s 2) m 2 ) X 1/2 (s 2 , ml M 2 ) . (46) 

ZS 2 

The limits on s 2 and t\ are given by the Chew-Low plot: for values of t\ bounded by 

tf = m 2 a + m\ - — [s + m 2 a - m 2 ] [s + m\- (m 2 + M) 2 ] 

± y/ /2 (s, m 2 a , m 2 b ) A 1/2 (s, (m 2 +M) 2 , m 2 ) , (47) 



the limits on s 2 are 



1 2 r 1 2 2l r 2 1 2 j. 1 



+ 



■^A^m^mgJA^.^X), (48) 



2m2 

s - = (m 2 + Mf . (49) 

Although not directly relevant for pions, we include for completeness the case where the 
hadronic masses m a and mi are unequal. In this case, t\ may become larger than given 
above, and the corresponding additional phase space is delimited by 

tf<h< (m a - mi ) 2 (50) 



and for this additional domain, s 2 is limited by s 2 of Eq. (|4q) and the corresponding lower 
branch, given by a minus sign before the A-factors. Clearly, t\ resembles the 2-phase space 
invariant t, but they are distinct quantities due to the presence of the additional momentum 

?■ 

Integration over pure phase space yields 

dR 3( ft X^js^mlM 2 ) 

[s,ti) = ds 2 = F(sJ) - F(s 2 ), (51) 

where 



dti Js 2 s 2 



F{s 2 ) = l 2 - {m\ + M 2 ) ln(s 2 + l 2 - m\ - M 2 ) 
- [m\ - M 2 ) In(s 2 + l 2 + m\ - M 2 ) 

-(M 2 -m 2 )ln(s 2 + / 2 -m 2 + M 2 ), (52) 

and we have written l 2 = A 1//2 (s 2 , m 2 , M 2 ) for short. This yields identical results to the 
pure phase space formula (|66D below, as it should. We note in passing that the ratio of the 
kinematic domains of (3-space) t\ to (2-space) t is exactly the phase space reduction factor 
found in (1691) below. 
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V. ANGLE-AVERAGED-CURRENT APPROXIMATIONS 



All the elements for our calculation are now in place. In this section, we look at three 
approximations of increasing sophistication. We proceed first to find explicit expressions 
for the pion-pion-photon cross section da^* of Eq. ( |25|) based on the factorization into 
a current and hadronic matrix element of Section [III A| . These approximations are based 
on an angular average of (— J 2 ) in the ah center-of-momentum frame. Using three-space 
invariants, we integrate the current covariantly in Section |V C. 



A. The Riickl formula 



In 1976, Riickl proposed j32| that the dilepton cross section be written in terms of the 
real photon cross section: 

rp p da h hi+e- en da hhl 



This ansatz, now in common use, was utilized by Haglin et al. p2j and others previously 
i7|,p!8|pT|,|3^| to derive dilepton production rates. In our notation, the Riickl formula reads 



a dr + dr_ 

da hhe+i - = 7r 2 M 2 dTq d(Thh ^ ' ( 54 ) 

which becomes, using the formulas of the previous section, and approximating / dVt + = Att 
and k ~ a/37r, 

do~hhe+e- = g K( ^ a Mvy • (55) 

The real-photon cross section da^ can be written down immediately from the virtual- 
photon cross section dahh-y* of (|25|), 

da hhl = 4™ (- J 2 )\M h \ 2 . (56) 

Neglecting q in the delta function of dR 3 , the 3-phase space integral of Eq. fl26|) factorizes 
into the usual 2-phase space and the qr-integration 



with 



leading all in all to 



dR 3 ~ dR 2 dr q , (57) 



dR 2 = 8(p a + p b -pi -p 2 ) dndr 2 , (58) 



da hh£+\~ 3 k a f 2 dr q 
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where 



da, 



iii, 



\M, 



, dR 2 
(2tt) 2 F 



(60) 



is the elastic hadronic cross section. Integrating (—J 2 ) first with respect to the photon phase 
space with dr q = dQ q \q\dq /2 leads to 



dr 

j2\ u,l Q 

' 2tt 



J*dq oy /<$-M*J 



dVL q 

Air 



-J 2 



(61) 



where A = [s + M 2 — (m± + m 2 ) 2 ]/2\fs is the upper kinematic boundary. 

The angular average over dQ q was found by Haglin et al. for the real-photon current ([16] 



in terms of velocities A in the (p a + p b = 0) cms as ^2 



Arc 



4 



-{Ql + QI + QI + Ql) - 2QaQbF ab - 2Q 1 Q 2 ^ 12 + 2^5] QiQjFij 

i=a,b j=l,2 



(62) 



where 



T 



2A 



In 



A-(/3j - A) - A. 



A-(A-A) + A. 



+ (*<-> i) 



(A -A) 2 - (Ax A 



1/2 



(63) 
(64) 



To convert these angular averages into functions of 2-phase space invariants (s,t), it is 
necessary to invoke approximation (C) restricting phase space: While p a + Pb = in the 
chosen frame, the presence of q means that pi + p 2 7^ 0, i.e. (3% and f3 2 are not true cms 
velocities. Ignoring therefore q at this point, the velocities and hence the angular averages 
can be written JEJ in terms of s, t and q , 



-I 2 ) 



4 



{Ql + Ql + QI + Ql) 

s — 2m 2 



-2(Q a Q b + Q 1 Q 2 
+ 2(Q a Q 1 + Q b Q 2 
+ 2(Q a Q 2 + Q 6 Qx 



s(s — Am 2 ) 
2m 2 -t 



In 



s + J s(s — 4m 2 ) 



t(s + u) 

2m 2 — u 
-u(s + t) 



In 



s — Js(s — Am 2 ) 
t(s + u) -t 



-t(s + u) +t 



In 



— u(s + t) — u 



—u(s + t) + u 



(65) 



with u = Am 2 — s — t as usual. 

Having expediently eliminated q at two points, it is very important in this formulation to 
correct for phase space, i.e. to alleviate the effects of approximation (C). Pure phase space 
integration of dR 3 (including q in the delta function) yields exactly 
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7TA 1 /2( s ' ;m 2 



m 



2\ 



dR 3 = J dr q v ^ b> (66) 

where 

s' = s - 2q oy fs + M 2 . (67) 
Integration over pure two-particle phase space, on the other hand, yields 

dR 2 =* Xl '\ mlml \ (68) 

so / dR,2, unconstrained by q, grossly overestimates / dR^, making it necessary to insert a 
posteriori the factor |36|] 

g(»,*) = ',Cy , i ,m, 'l < 69 > 

s' X i ' 2 (s, mf, m^) 

into the go-integration fl6l|) . The "Riickl approximation" for dilepton production therefore 
reads, in its final form, 

= \^~J d ^(s, t) J dq \q\ C(s, g ) <-/ 2 (*, t)> ; (70) 



this is the form used in Ref. [221. 



B. Current formulae based on Lichard 



Starting not from the Riickl formula but rather from Lichard's improved formalism |26 
of Section [111 A| , the derivation proceeds along the same lines as the above. One finds 
correspondingly for the real-photon current 



drr L0 

au hM+l 



) k a 
dM 2 ' ~ JPn 



J da hh {s, t) J dq \q\ C{s, q ) (~I 2 {s, t)) , (71) 



differing from the above by a constant factor 3/2 and the extra factors entering the exact 
version (|24[) of k. 

The corresponding virtual-photon current yields the third approximation within the 
angle-averaged current family, 

d ~%f^ = jjS- J 'da hh (s,t) Jdq \ q \C(s,q )(-J 2 (s,t)). (72) 
where the angular average of (pT5|) is 

(Ql + Q 2 b )(4m 2 -M 2 ) 



qls - q 2 (s - Am 2 ) - 2q ^M 2 + M 4 
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where 



qls — q 2 (& 



+ 



+ 



(Q 2 + Q 2 )(Am 2 -M 2 ) 



Am 2 ) + 2q^sM 2 + M 4 



_ Q a Q b {2s - Am 2 - M 2 ) ^ |g|Vg-4m 2 + g 0v ^-M 2 

\q\Vs-Am 2 -q y/s + M 2 
\q\y/s- Am 2 + q ^/s + M 2 



\q\v-ys — Am 2 
Q 1 Q 2 (2s - Am 2 - M s 

'— i -i 1 - : | : 

, \q\y/s - Am 2 - q ov fs - M 2 

{Q a Q 1 + Q b Q 2 ){Am 2 -2t + M 2 ), 
m 



-In 



(73) 



2 X t 

(QaQ 2 + Q fc Qi)(4m 2 - 2u + M 2 ) 
2Xu 



In 



q 2 t - X t- M 2 v_)(q 2 t - X t + M 2 v+) 
q 2 t + X t- M 2 v^){q 2 t + X t + M 2 v + ) 
q 2 u -Xu~ M 2 v_)(q 2 u - X u + M 2 v + ) 
! ' X u - M 2 v.)(q 2 u + X u + M 2 v + ) 



(q 2 u 



v± = g v^ ± M 2 , 

Xt = \q\^ -t{sq 2 + uq 2 ) - uM± 

Xu = \q\J-u{sq 2 + tq 2 )-tM*. 



(74) 
(75) 
(76) 



Eq. (|73|) reduces to Eq. ( |65| ) for M — > 0. We shall be testing the effect of using the virtual 



over real-photon currents within the present approximation in Section [VIII| below. We note 
in passing that Eq. (|73| ) would be preferable to Eq. (|65D for very small hadron masses m, 
since then the latter is logarithmically sensitive to 1/m while the former is regulated by the 
dilepton invariant mass M. 



C. Three-phase space current approximations 

Three approximations were needed for the above angle-averaged current cross sections: 
the Riickl formula, converting cms velocities f3i to invariants s and t leading to Eqs. (|62D and 
(IH), and correcting for phase space using Eq. (|69|) . At the cost of increased computing time, 
all three can be avoided by utilizing the formalism of three-particle phase space invariants. 

Within 3-phase space, and writing the flux as F = 2X 1 ^ 2 (s,m 2 ,m 2 ) : the bremsstrahlung 
cross section of Eq. fl25]) becomes[] 

duhhi+i- Ana k(M 2 )7t r , . 2 f ds2ds\dt2 



dM 2 w (2vr) 5 M 2 8A(s, m 2 ,m 2 ) J 11 ftV ' i;i J s/b 



s ) ~ /o \**r> oT7 5 2\ dh\M h (s,ti)\ / 7= —J (s,ti,s 2 , si,t 



'2, 



(77) 

This differs from Eq. (|72|) in the way the current is treated: by an angular average plus 
go-integration for the earlier, by full 3-phase space invariants in the present case. Note that 
M.h is here a function of t\ rather than t. 



1 Technically, it would in equation ([77]) be possible to restore q to Aih and then to include it 
with — J 2 in the inner integral, but because of missing contact terms this procedure is not gauge 
invariant and shall not be pursued. 
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VI. EXACT CROSS SECTION 



In all the above approximations, neglecting the dependence on q of the squared hadronic 
matrix element \Mh(Pa,Pb,Pi,P2)\ 2 is unavoidable because this had been a precondition to 
factorizing into electromagnetic and hadronic parts in Eq. (H). There have been attempts 
go beyond this on-shell matrix element approximation by expanding hA h in q (see Refs. 
26| , |27|j3^| ) , but this procedure involves successively more complicated derivatives without 



guaranteeing convergence. 

Using 3-phase space invariants and the exactf] expressions ( j36[) for the pion-pion-photon 
cross section and ( |5oD for the dilepton cross section, such expansions in q become superfluous. 
To calculate explicit dilepton cross sections, we merely insert from our favorite model for 
pion-pion scattering the complete set of matrix elements, and write the elements of the 
squared sum of fl36|) in terms of these invariants: 



M^M* niM (p a ,Pb,PuP2,q) — ► M'^M* ni ,(s,t 1 ,s 2 ,s 1 ,t 2 ) =T mn . (78) 

This casting in terms of invariants means that q and all consequent off-shell effects in \M.\ 2 
are fully taken care of] 3 ] 



Putting it all together, we obtain 
d a hhe+e- 47ra k,(M 2 )tc r dt 1 ds 2 dsidt 2 



dM 2 w (2tt) 5 M 2 8\(s,m 2 a ,m 2 ) J 



-^T mn (s,ii,S2,si,t 2 ) 



inn 



(79) 



While the cross sections may be exact, microscopic models for the pion-pion interaction are 
far from perfect. In the next section, we make use of a simple model to illustrate the present 
cross section and the approximations enumerated before. The main point of this section, 
however, is to stress that, apart from model-induced ailments, the exact treatment of any 
model for the reaction titt — > 7ttt£ + £~ is possible. 



VII. GAUGE-INVARIANT OBE MODEL 

To compare dilepton cross sections under the exact formula and various approximations 
introduced so far, we must make a model: little if any experimental data for the differential 
pion-pion-photon cross section is available. To this purpose, we return to a simple lagrangian 



with a, p and /(1270) exchange used previously [22| which, when t-integrated, successfully 



reproduced experimental data for elastic tx + h — > n + n scattering [ ft7| , |35| ]: 

A„t = g a <rd^-d»iT + gp^TzKd^ + g f f^ l n-d v iz . (80) 



2 To order a for the emission of real photons and to order a 2 for dilepton emission. 

3 Matrix elements keeping their full dependence on q can be implemented within 2-phase space 
when only the differential cross section da/d 3 qdM 2 need be calculated; see e.g. Ref. j29[] ; for 
da/dM 2 , this becomes impossible. 
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Here, we include the full momentum dependence of the / propagator which now reads 



39 



■ 40fl 



aS9/3-/ 9a 

kfikd kpkj k a k$ k a k^ \ 

mj mj mj mj J 

1 k a kp\ (l k 7 k s \ 




The p field also has a fc-dependent term in the unitary gauge propagator |41 
iR^ v (k) 



k 2 — m 2 + im T, 



kuky 

-9iw + 



m 2 



(83) 



p J 

but this term contributes to matrix elements only in the form (m 2 — m\)(m 2 — m 2 ,) etc. and 
is thus zero for equal pion masses. 

For t- and w-channel exchange of an a-meson (a = a, p /) with momentum k by the 
pions, we implement monopole strong form factors 

M*) = (84) 
In Figure 4, we show the results of fitting this model to the elastic scattering data as before 



1 22|| but with the different propagator (|82|). Again, the fit is reasonably good, but the 
strong dependence on momentum-dependent terms in the / propagator shows up in the 
region above the / peak. Best-fit values for the free parameters now are: m a = 475 MeV, 
m p = 775 MeV, m/ = 1220 MeV; m a g a = 3.1, g p = 6.15, m/<7/ = 8.85. Throughout, we 
take m n = 140 MeV. 

For the current-based approximations of Section [V|, the above lagrangian is all we need; 
for the complete cross section of Section [VI|, it must be augmented by internal photon 
emission pieces. This can be achieved by minimal substitution, thereby guaranteeing gauge 
invariance (for a similar approach, see the "Quantum Hadrodynamics" lagrangians of Ref. 
T[). To the hadronic lagrangian (RO) we therefore add 

Details regarding these pieces and corresponding vertex factors can be found in the Ap- 
pendix. 

The ppj interaction results in a modified strong form factor for t- or u-channel charged 
p exchange. If h k and hi are the usual form factors at the 7mp(k) and irirp(l) vertices 
respectively, with k = I + q, electromagnetic gauge invariance requires that the form factor 
for the emission of a photon by the p be 

H(kj)= ^- Vh l , (86) 
[u - v)h k hi 
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where u = k 2 — m 2 + im p Y p and v = I 



the monopole form (|84f) of Ai this becomes j42l , [43 



rri^+ im p Tp are the respective p propagators. For 



ff(M) = 1 + 



A- 2 



im p Y p s 



P 



m. 



I 2 — im n T 



A' 2 



(87) 



Note that we do not include strict Vector Meson Dominance coupling of the photons to 
vector mesons. This is equivalent to the introduction of electromagnetic form factors. For 
the specific application we are considering, small M's, those would be ~ 1. For previous 



OBE models, see e.g. Refs. [2S , 44] and references therein 



VIII. RESULTS AND CONCLUSIONS 

To summarize: the differential cross section for dilepton production as a function of 
invariant mass, da^hi+e- /dM, can be found in a number of approximations and an exact 
way (within the limitations of the hadronic reaction model). The approximations are: the 
Ruckl approximation ([70| ), two approximations, Eqs. (|7T|) and (|72|), based on a more care- 
ful derivation by Lichard, using the real and virtual photon currents respectively, and an 
approximation integrating the current covariantly in 3-phase space, Eq. ([FT]). The "exact" 
formula is given by Eq. (|79|) . 

Of these, the Ruckl and Lichard cross sections can be found easily and quickly within 
2-phase space. The 3-phase space current approximation is somewhat more difficult but 
still uses 2-phase space hadronic matrix elements. Implementing the exact formulation ( |79|) 
involves much more effort: if there are N diagrams contributing to a given reaction, the sum 
in ([75|) contains N(N + l)/2 terms. Since an individual term T mn in ( [78] ) is itself a Dirac 
sum, only partial factorization of the overall sum is possible. Further complications arise 
from the presence of imaginary pieces in the matrix elements. 

To quantify the differences between the approximations and the exact formulation, we 
have studied all five distinct pion-pion reactions. Writing (H — ) — > (H — ) as shorthand for 
the reaction (tt + tt~ — > tt + tt~£ + £~) and so on, we have calculated within our OBE model cross 
sections for (+-) -> (+-), (++) -> (++), (+-) -> (00), (00) -> (+-) and (+0) -> (+0). 
A total of 36 diagrams contribute at tree level to the first two reactions (24 for 7* emitted 
by external pion lines, 12 by vertices and the exchange meson); the remaining reactions are 
made up of 16 diagrams each (8 external, 8 internal). 

Numerical results were checked by performing the following consistency checks: setting 
Si — > s, t\ — > t, t 2 — > ml, s 2 — > ml and M 2 — ► in the hadronic part of the exact matrix 
elements (while keeping these variables where they enter the "current" part) and using only 
the external-emission diagrams reproduces exactly the results of the 3-phase space integrated 



current (|77|) . Gauge invariance in the a, p and / fields provided another sensitive test. 

Figures 5-9 show the cross sections for these five reactions as a function of ^/s. In 
every case, the left and right panels show da/dM for M = 10 MeV and 300 MeV respec- 
tively. Final-state symmetrization factors were included where appropriate. Initial-state 
symmetrization was also included for (++) — > (++) and (00) — > (H — ) in order to facilitate 
use within a thermal pion gas environment. When the initial-state pions are considered dis- 
tinguishable, these two cross sections should be multiplied by a factor 2 for these reactions. 
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Because they are identical in structure to their charge-conjugate versions, cross sections for 
the reactions (+0) — > (+0) and (++) — > (++) were doubled — this, too, should be corrected 
for when desired. 

To prevent overcrowding, only the Ruckl plus the virtual-7* current approximations (|72|) 
and ( |TTD are shown, together with the exact calculations. All these were computed using 



the same OBE model and parameter values specified in Section [VI I| above. 

In the complicated structures and deviations in Figs. 5-9, the following points are of 
interest : 

1. The p peak in the reaction (H — ) — ► (H — ) is overestimated by a factor 2.5 and 1.5 by the 
Ruckl and Lichard approximations respectively for M — 10 MeV; for M = 300 MeV, 
the overestimation is much greater (Figure 5). Only the 3-space current approximation 
does an adequate job here. Similarly, the / peak of the exact result lies well below the 
corresponding approximations. Overestimation (1.7-5.0) of the p peak also occurs for 
approximations in (+0) — > (+0). 

2. In the reaction (+0) — > (+0), the 3-space current approximation, on the other hand, 
lies well below the exact result (Figure 6). 

3. Overestimation factors for da/dM are of order 0.5-4 for the reactions (00) — > (H — ) 
and (-1 — ) — > (00), depending on approximation and cms energy. Here, too, the 3- 
space current approximation generally does better than the others but again tends to 
underestimate in some parts (Figs. 7, 8). 

4. The largest discrepancy between approximations and the exact result occur for the 
reaction (++) — > (++) (Figure 9): for the Ruckl approximation, factors 3 (for 10 
MeV) to 30 (for 300 MeV) arise, while the Lichard approximation yields corresponding 
overestimation factors of 1.9 and 14-20. Since, however, this reaction contributes only 
little when contributions from all reactions are added up, this effect is not of much 
interest in a heavy-ion context. 

5. There is generally not much difference between using the real-photon (not shown) or 
virtual-photon currents in the Lichard approximations; this is in agreement with the 



conclusions of Ref. p7j- For M = 300 MeV, some differences can be discerned, but 



both deviate in general more from the exact result than from each other. 

Figures 10 and 11 show the dilepton production rates per unit spacetime summed over all 
seven reactions, calculated using the Boltzmann formula fl3"7p, for temperatures T = 100 and 
200 MeV respectively. These temperatures probably lie below and above an expected transi- 
tion to a new quark-gluon phase and so represent extrema in their behavior. Corresponding 
ratios of approximation over exact rates are shown in Figures 12 and 13. The Ruckl approx- 
imation is the worst, as expected; overestimation is in the range 2-4 for T = 100 MeV, and 
2-8 for T = 200 MeV. The Lichard approximations overestimate by factors 1.4-4, depending 
on temperature and s. Here, the real-'-) 3-phase space current, the lower of the two dotted 
lines, does surprisingly well, deviating from the exact result by less than 20% throughout. 
Comparing this fact, however, with the substantial deviations shown in the s-differential 
cross sections of Figs. 5-9, one must conclude that much of the "agreement" seen on the 
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level of rates is due to averaging over differences seen in the s-dependent cross sections and 
that such "agreement" therefore does not validate the approximations. 

We also note the fact that none of the approximations approaches the exact result for 
small values of M: even for the smallest value shown (M = 10 MeV), the discrepancy is 
still above 40% for the Lichard approximations and larger than a factor 2 for the Riickl 
approximation. 

Figure 14 shows the fractional contributions to the total T = 200 MeV rate for the 
Riickl, virtual-7 current Lichard, and exact cross sections. As discussed, the (++) — > (++) 
contribution virtually falls away in the exact calculation; by contrast the reaction (+0) — > 
(+0) plus charge conjugate becomes more important, in second place behind the dominant 
reaction (H — ) — > (H — ). 

Finally, we ask where the discrepancies between approximations and exact cross section 
originate: is this due to the inclusion of diagrams with photon emission from vertices and 
propagators? A partial answer is provided by Figure 15, where the cross sections for reactions 
(++) — > (++) and (H — ) — > (00) are plotted. The upper lines correspond, as before, to the 
Riickl and two Lichard approximations, while the solid line again represents the exact result. 
The lowest dash-dotted line, on the other hand, represents the exact result but excluding 
all internal diagrams and their cross terms with external ones. The difference between 
this lower line and the exact result (solid line) therefore represents the contribution of the 
internal diagrams; while the difference between the lower dash-dotted line and the upper 
lines (approximations) represents the change in cross section due to inclusion/exclusion of 
q in the external-emission diagrams. 

We see that the contribution of internal emission is not all that large, albeit nonnegligible. 
By far the most important effect on da/dM is the inclusion of the full dependence of Ai on 
the photon momentum q. In fact, the effect is so large that for the reaction (+- 1-) — > (++) 
the cross section would become negativefj if internal diagrams were not included! In other 
words, among the list of approximations listed in Section [11 Bj , approximation (B) is the 
most far-reaching. 

The importance of including q fully in the matrix element can be understood quali- 
tatively: as illustrated in Figs. 5-9 for one reaction, the elastic pion-pion cross sections 
vary considerably within the region of interest. Adding or subtracting a sizeable amount 
of energy-momentum in the form of q therefore naturally leads to a substantial change in 
cross section also. This change is, however, neglected in all the approximations listed^. We 
believe that this (unavoidable) neglect is at the heart of the considerable differences between 
approximations and exact result seen in Figs. 5-9. 

In summary, the various approximations for bremsstrahlung from pion-pion collisions 



Since this curve does not represent a gauge- invariant calculation, this negative cross section is 
of no physical consequence. In fact, it illustrates once again the necessity of including internal 
diagrams to make the cross section positive again. 

5 This also explains why the overestimation is so large for the reaction (++) — ► (++): there are 
no resonances in the cross section which could provide some amount of cancellation by pion-pion 
reactions being "shifted into" and "shifted out of" a resonance peak. 
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cannot in general be believed beyond at best a factor 2 or more, depending on invariant 
mass and temperature. 

We have shown in Sections [I I [13] and [VI] how 3-phase space can be fully taken into 
account. The most important effect, it turns out, is not the inclusion of photon emission 
from vertices or propagators but the strong dependence of the cross section on the inclusion 
of the photon momentum q into the hadronic collision part. 

Since this calculation neglected baryonic degrees of freedom and resonances, and because 
of the technical difficulties of a further integration over spacetime, our results are not applied 
directly to experimental data; this is a matter for the future. All we can say at this stage 
is that it is likely that the approximations overestimate the dilepton yield throughout the 
mass region of interest, and that the degree of overestimation rises both with increasing 
invariant mass M and temperature T. Since the bremsstrahlung channel competes with 
various Dalitz decays, and since various experiments take much care in applying cuts to 
eliminate Dalitz pairs as far as possible a quantitative estimate would necessarily have 
to take such experimental cuts into account also. 

In addition, attention should be paid to the issue of the Landau- Pomeranchuk effect, 
which is expected to suppress dilepton bremsstrahlung rates p^j25|j28[1 . Our results, taken 



in conjunction with these calculations, therefore appear to indicate even larger suppression 
of bremsstrahlung dileptons than had been thought previously. 
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APPENDIX: FULL QBE LAGRANGIAN 



Besides the hadronic interaction lagrangian fl80|) with its vertex factors 



wrrp 



7T7T/ 



2ig a p ■ p' 

ig* P ■ p' 
Tig P {p + p'Y 

±ig P (p + p'Y 
±ig p {p + p'Y 
igf(p a p ,f3 + pfip la ^ 
ig f p a p'P 



for TT^ij)) 

for 7T°(p) ■ 
for ^(p) 
for ^(p) 
for 7i° (p) ■ 
for TT^ijp) 
for 7f°(p) ■ 



' 7f ± {p') cr 
ttV) o 

•*\p')pI 

7T°{p') f a ? 



we need the electromagnetic interaction lagrangian linear in A^, 
With the pion current 



(Al) 



(A2) 
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# = i(7r-^7T + - TT+^TT") , (A3) 

the p field tensor, 

and F pu = <9 M Aj, — (9^A M the usual electromagnetic tensor, the individual contributions are 



r - —pA 



7T7TCr7 



TVTVpJ 



ieA, [ P - v {dP P v + - d v p\) - p +v {d»p»_ - d v pt) 
-2eg a A^ , 
-egpAp [7rx(7rxp^)] 3 
eg p A^ 2ty + 'k~Pq — 7r~7T°p+ — 7r + 7r°p- 
-egffr Wj: + A" ft] . 



ied ^ A v 



U V V U 

P-P+ - P-P+ 



£-/r 71-/7 

The corresponding vertex factors are 



(A5) 
(A6) 

(A7) 
(A8) 

(A9) 



7T7r7 



PfH 



wniry 



wnp'y 



7T7r/7 



g af5 (p + p'Y 
±2ieg a (p + p'Y 
-2ieg p g^ 
+ieg p g> M/ 
+ieg p g^ u 
±iegf g afl (p + p'Y + g^ijp + p') 



for n ± (p) 
for p±{p) 
for 7r ± (^?) 
for ^(p) 
for ^(p) 



pl(p')A»(q) 

7T ± (p')^ 
7T°(p')P±^ 



(A10) 



for 7T°(p)p^ -> 7T ± (p')^ 
for 7T ± (p) -> 7T ± (p')/ a/3 ^ 
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List of Figures 

Figure 1: Breakdown of the soft photon approximation: the SPA is valid only when q 
is much smaller than the two solid lines shown. The region between the dashed lines is 
the domain of integration for q when calculating cross sections as a function of dilepton 
invariant mass. 

Figure 2: One of the four contributing diagrams following from SPA factorization of the 
matrix element into leptonic tensor L^, photon emission current J M , and elastic hadronic 
matrix element Aih- 

Figure 3: Relativistic invariants for three-particle final state phase space. 

Figure 4: OBE model fitted to the elastic 7r + 7r~ — > 7r + 7r~ cross section. Parameter values 
determined from this fit are used for all reactions throughout this paper. 

Figure 5: Cross section for the reaction 7r + 7r~ — > 7r + 7r~e + e~ as a function of yfs for fixed 
dielectron invariant masses M = 10 and 300 MeV. Solid line: exact OBE calculation Eq. 
([79]). Dash-dotted line: Riickl approximation (|70|) — this is the commonly-used version of 
the SPA. Other approximations are: Dashed line: angular-averaged virtual 7 current (f?"2"|). 
Dotted line: 3-phase space current (|77|). 

Figure 6: Same as Figure 5, for the reactions 7r + 7r° — > n + iT e + e~ plus 7r~7r° — > n^n°e + e~ 

Figure 7: Same as Figure 5, for the reaction 7r + 7r~ — > 7r°7r°e + e~, including final-state sym- 
metry factor. Note the different scales on the y-axes. 
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Figure 8: Same as Figure 5, for the reaction 7r°7r — > 7r + 7r e + e . An initial-phase space 
factor of 1/2 is included. 

Figure 9: Same as Figure 5, for the reactions 7r + 7r + — > 7r + 7r + e + e~ plus n~7T~ — > 7r~7r~e + e~. 
Initial- times final-state symmetry factors (1/2) 2 are included. 

Figure 10: Total bremsstrahlung yield of e + e~ for all seven pion-pion reactions, for a Boltz- 
mann gas with temperature T = 100 MeV. Solid line: exact OBE calculation, Dash-dotted 
line: Riickl approximation, Dashed line: angular-averaged virtual-7* current, Dotted line: 
3-phase space current, virtual-7* current (upper) real-7 current (lower). 

Figure 11: Same as Figure 10, for temperature T = 200 MeV. 

Figure 12: Ratios of SPA approximation calculations divided by exact OBE rate for 
T = 100 MeV. Dash-dotted line: Riickl/exact, Dashed line: (angular-averaged virtual 7 
current) /exact. Upper dotted line: (3-phase space virtual 7 current) /exact. Lower dotted 
line: (3-phase space real 7 current) /exact. 

Figure 13: Same as Figure 12, for T = 200 MeV, 

Figure 14: Fractional contribution to the total dielectron rate at T = 200 MeV as a function 
of dilepton invariant mass M for (a) the Riickl, (b) virtual-photon Lichard approximations 
and (c) the exact calculation. Solid line (upper): (H — ) — > (H — ), Dashed line: (+0) — > (+0) 
plus c.c, Solid line (lower): (00) -> (H — ), Dash-dotted line: (H — ) -> (00), Dotted line: 
(++) -> (++) plus c.c. 

Figure 15: Contribution of external-emission vs. external-plus-internal diagrams for the 
reactions (+-) -> (00) (left) and (++) -> (++) (right), both for M = 300 MeV. Lines 
are as in Figs. 7 and 9. The new dash-dotted line below the (solid line) exact calculation 
represents contributions arising solely from emission of 7* by an external pion line, but 
taking q into account in A4, in contrast to the approximations (upper lines) which neglected 
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